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Magneto-transport in mesoscopic rings and cylinders: Effects of electron-electron 

interaction and spin-orbit coupling 

Santanu K. Maiti^'B 
'School of Chemistry, Tel Aviv University, Ramat-Aviv, Tel Aviv-69978, Israel 

We undertake an in-depth analysis of the magneto-transport properties in mesoscopic single- 
channel rings and multi-channel cylinders within a tight-binding formalism. The main focus of this 
review is to illustrate how the long standing anomalies between the calculated and measured current 
amplitudes carried by a small conducting ring upon the application of a magnetic flux (p can be 
removed. We discuss two different cases. First, we examine the combined effect of second-neighbor 
hopping integral and Hubbard correlation on the enhancement of persistent current in presence of 
disorder. A significant change in current amplitude is observed compared to the traditional nearest- 
neighbor hopping model and the current amplitude becomes quite comparable to experimental 
realizations. In the other case we verify that in presence of spin-orbit interaction a considerable 
enhancement of persistent current amplitude takes place, and the current amplitude in a disordered 
ring becomes almost comparable to that of an ordered one. In addition to these, we also present the 
detailed band structures and some other related issues to get a complete picture of the phenomena 
at the microscopic level. 

PACS numbers: 73.23.-b, 73.23.Ra, 73.21.Hb 



I. INTRODUCTION 

The study of magneto-transport properties in low- 
dimensional systems is always an interesting problem 
in mesoscopic physics. This is relatively a new branch 
of condensed matter physics which deals with systems 
whose dimensions are intermediate between the micro- 
scopic and macroscopic length scales^"— . In contrast to 
the macroscopic objects where we generally use the laws 
of classical mechanics, the meso-scale systems are treated 
quantum mechanically since in this region fluctuations 
play the very crucial role. Over the last many years low- 
dimensional model quantum systems have been the ob- 
jects of intense research, both in theory and in experi- 
ments, mainly due to the fact that these simple looking 
systems are prospective candidates for nano devices in 
electronic as well as spintronic engineering^^—. 

Several striking spectral properties are also exhibited 
by such systems owing to the quantum interference which 
is specially observed in quantum geometries with closed 
loop structures. The existence of dissipationless current 
in a mesoscopic metallic ring threaded by an Aharonov- 
Bohm (AB) flux (/> is a direct consequence of quantum 
phase coherence. In this new quantum regime, two im- 
portant aspects appear at low temperatures. The first 
one is that the phase coherence length L^ i.e., the length 
scale over which an electron maintains its phase memory, 
increases significantly with the lowering of temperature 
and becomes comparable to the system size L. The other 
one is that the energy levels of such small finite size sys- 
tems are discrete. These two are the most essential cri- 
teria for the appearance of persistent charge current in 
a small metallic ring/cylinder due to the application of 
an external magnetic flux 0. In the pioneering work of 
Biittiker, Imry and Landauei— , the appearance of persis- 
tent current in metallic rings has been explored. Later, 
many excellent experiments^^"— have been carried out 



in several ring and cylindrical geometries to reveal the 
actual mechanisms of persistent current. Though much 
efforts have been paid to study persistent current both 
theoretically^"— as well as experimentally^^"— , yet sev- 
eral anomalies still exist between the theory and experi- 
ment, and the full knowledge about it in this scale is not 
well established even today. 

The results of the single loop experiments are signifi- 
cantly different from those for the ensemble of isolated 
loops. Persistent currents with expected 0o periodic- 
ity have been observed in isolated single Au rings^^ and 
in a GaAs-AlGaAs ring21. Levy et al^ found oscilla- 
tions with period 0o/2 rather than 0o in an ensemble 
of 10^ independent Cu rings. Similar 4>o/2 oscillations 
were also reported for an ensemble of disconnected 10^ 
Ag rings2S as well as for an array of 10^ isolated GaAs- 
AlGaAs rings^. In a recent experiment, Jariwala et al^ 
obtained both (f>o and (/>o/2 periodic persistent currents in 
an array of thirty diffusive mesoscopic Au rings. Except 
for the case of the nearly ballistic GaAs-AlGaAs ring22, 
all the measured currents are in general one or two orders 
of magnitude larger than those expected from the theory. 

Free electron theory predicts that at absolute zero tem- 
perature (T = OK), an ordered one-dimensional (ID) 
metallic ring threaded by magnetic flux (j) supports per- 
sistent current with maximum amplitude /q = evp/L, 
where vp is the Fermi velocity and L is the circumfer- 
ence of the ring. Metals are intrinsically disordered which 
tends to decrease the persistent current, and the calcu- 
lations show that the disorder-averaged current (/) cru- 
cially depends on the choice of the ensemble^i^i^. The 
magnitude of the current (J^)^" is however insensitive to 
the averaging issues, and is of the order of IqI/L, I being 
the elastic mean free path of the electrons. This expres- 
sion remains valid even if one takes into account the finite 
width of the ring by adding contributions from the trans- 
verse channels, since disorder leads to a compensation 



between the channels^i^. However, the measurements 
on an ensemble of 10^ Cu rings^^ reported a diamagnetic 
persistent current of average amplitude 3 x 10"'^ evpjL 
with half a flux-quantum periodicity. Such 0o/2 oscilla- 
tions with diamagnetic response were also found in other 
persistent current experiments consisting of ensemble of 
isolated rings^Si^. 

Measurements on single isolated mesoscopic rings on 
the other hand detected (/>o-periodic persistent currents 
with amplitudes of the order of /p ^ evp / L, (closed to 
the value for an ordered ring). Theory and experiment^ 
seem to agree only when disorder is weak. In another 
recent nice experiment Bluhm et al^ have measured 
the magnetic response of 33 individual cold mesoscopic 
gold rings, one ring at a time, using a scanning SQUID 
technique. They have measured h/e component and pre- 
dicted that the measured current amplitude agrees quite 
well with theory^ in a single ballistic ring^i and an en- 
semble of 16 nearly ballistic rings^. However, the am- 
plitudes of the currents in single-isolated-diffusive gold 
rings2^ were two orders of magnitude larger than the 
theoretical estimates. This discrepancy initiated intense 
theoretical activity, and it is generally believed that the 
electron-electron correlation plays an important role in 
the disordered diffusive rings^"— . An explanation based 
on the perturbative calculation in presence of interaction 
and disorder has been proposed and it seems to give a 
quantitative estimate closer to the experimental results, 
but still it is less than the measured currents by an order 
of magnitude, and the interaction parameter used in the 
theory is not well understood physically. 

To remove the controversies regarding the persistent 
current amplitude between theoretical and experimen- 
tal verifications we can proceed in two different ways. 
In literature almost all the theoretical results have been 
done based on a tight-binding (TB) framework within the 
nearest-neighbor hopping (NNH) approximation. It has 
been shown that in the NNH model electronic correla- 
tion provides a small enhancement of current amplitude 
in disordered materials i.e., a weak delocalizing effect is 
observed in presence of electron-electron (c-c) interac- 
tion. As a first attempt, we modify the traditional NNH 
model by incorporating the effects of higher order hop- 
ping integrals, at least second- neighbor hopping (SNH), 
in addition to the NNH integral. It is also quite physical 
since electrons have some finite probabilities to hop from 
one site to other sites apart from nearest-neighbor with 
reduced strengths. We will show that the inclusion of 
higher order hopping integrals gives significant enhance- 
ment of current amplitude and it reaches quite closer to 
the current amplitude of ordered systems. This is one ap- 
proach. In the other way we examine that in presence of 
spin-orbit (SO) interaction a considerable enhancement 
of persistent current amplitude takes place, and the cur- 
rent amplitude in the disordered ring is almost compa- 
rable to that of an ordered ring. The spin-orbit fields 
in a solid arc called the Rashba spin-orbit interaction 
(RSOI) or the Drcsselhaus spin-orbit interaction (DSOI) 



depending on whether the electric field originates from 
a structural inversion asymmetry or the bulk inversion 
asymmetry respectively^. Quantum rings formed at the 
interface of two semiconducting materials are ideal can- 
didates where the interplay of the two kinds of SOI might 
be observed. A quantum ring in a heterojunction is real- 
ized when a two dimensional gas of electrons is trapped in 
a quantum well due to the hand offset at the interface of 
two different semiconducting materials. This hand offset 
creates an electric field which may be described by a po- 
tential gradient normal to the interface^. The potential 
at the interface is thus asymmetric, leading to the pres- 
ence of a RSOI. On the other hand, at such interfaces, 
the bulk inversion symmetry is naturally broken. 

The other important controversy comes for the deter- 
mination of the sign of low-field currents and still it is 
an unresolved issue between theoretical and experimen- 
tal results. In an experiment on persistent current Levy 
et al^ have shown diamagnetic nature for the measured 
currents at low- field limit. While, in other experiment 
Chandrasekhar et al^ have obtained paramagnetic re- 
sponse near zero field limit. Jariwala et al^ have pre- 
dicted diamagnetic persistent current in their experiment 
and similar diamagnetic response in the vicinity of zero 
field limit were also supported in an experiment done by 
Deblock^S et al. on Ag rings. Yu and Fowler— have 
shown both diamagnetic and paramagnetic responses in 
mesoscopic Hubbard rings. Though in a theoretical work 
Cheung et al^ have predicted that the direction of cur- 
rent is random depending on the total number of elec- 
trons in the system and the specific realization of the 
random potentials. Hence, prediction of the sign of low- 
field currents is still an open challenge and further studies 
on persistent current in mesoscopic systems are needed 
to remove the existing controversies. 

In the present review we address several important is- 
sues of magneto-transport in single-channel mesoscopic 
rings and multi-channel mesoscopic cylinders which are 
quite challenging from the standpoint of theoretical as 
well as experimental research. A brief outline of the pre- 
sentation is as follows. 

First, we address magnetic response in mesoscopic 
Hubbard rings threaded by AB flux (f>. We try to propose 
an idea to remove the unexpected discrepancy between 
the calculated and measured current amplitudes by in- 
corporating the effect of second-neighbor hopping (SNH) 
in addition to the traditional nearest-neighbor hopping 
(NNH) integral in the tight-binding Hamiltonian. Us- 
ing a generalized Hartree-Fock (HF) approximation^"—, 
we numerically compute persistent current (/), Drude 
weight (D) and low-field magnetic susceptibility (x) as 
functions of AB fiux 0, total number of electrons N^ and 
system size N. With this (HF) approach one can study 
magnetic response in a much larger system since here 
a many-body Hamiltonian is decoupled into two effec- 
tive one-body Hamiltonians. One is associated with up 
spin electrons and other is related to down spin electrons. 
But the point is that, the results calculated using gener- 



alizcd HF mean-field theory may deviate from exact re- 
sults with the reduction of dimensionality. So we should 
take care about the mean-field calculation, specially, in 
ID systems. To trust our predictions, here we also we 
make a comparative study between the results obtained 
from mean-field theory and exactly diagonalizing the full 
many-body Hamiltonian. The later approach where a 
complete many-body Hamiltonian is diagonalized to get 
energy eigenvalues is not suitable to study magnetic re- 
sponse in larger systems since the size of the matrices 
increases very sharply with the total number of up and 
down spin electrons. 

Next, we explore the behavior of persistent current in 
an interacting mesoscopic ring with finite width threaded 
by an Aharonov-Bohm flux 0. For this cylindrical sys- 
tem we also see that the inclusion of higher order hopping 
integrals leads to a possibility of getting enhanced per- 
sistent current and the current is quite comparable to 
the ordered one. Our results can be utilized to study 
magnetic response in any interacting mesoscopic system. 

Finally, in the last part, we focus our attention on the 
behavior of persistent current in a one-dimensional meso- 
scopic ring threaded by a magnetic flux in presence of 
the Rashba and Dresselhaus SO interactions. Here, the 
effect of electron-electron interaction is neglected. We 
show that the presence of the SO interaction leads to a 
significant enhancement of persistent currenlj^. In ad- 
dition to these, we also describe very briefly the energy 
band structures and the oscillations of persistent current 
as the RSOI is varied to make the present communication 
a self contained study. 

Throughout the review we perform all the essential fea- 
tures of magneto-transport at absolute zero temperature 
and set c = e = ft, = 1 for numerical calculations. 



II. A HUBBARD RING IN ABSENCE OF SO 
INTERACTIONS 

In this section we describe the magneto-transport 
properties in a single-channel ID mesoscopic ring in the 
presence of on-site Coulomb interaction. The effect of 
SO interaction is not taken into account. 



A. Model and theoretical formulation 

We start by referring to Fig. [1] where a normal metal 
ring is threaded by a magnetic flux ^. To describe the 
system we use a tight-binding framework. For a A^-site 
ring, penetrated by a magnetic flux (measured in unit 
of the elementary flux quantum 00 = ch/e)^ the tight- 
binding Hamiltonian in Wannier basis looks in the form, 

i.(T ij,a 

ik,a i 



where, eic, is the on-site energy of an electron at the site i 
of spin a (t, l). The variable t corresponds to the nearest- 
neighbor (j ~ i± 1) hopping strength, while ii gives the 
second-neighbor (k = i±2) hopping integral. 6 = 2tt(J)/N 
and 6*1 = An4)/N are the phase factors associated with the 
hopping of an electron from one site to its neighboring 



site and next-neighboring site, respectively. 



and 




FIG. 1: (Color online). Schematic view of a ID mesoscopic 
ring penetrated by a magnetic flux (f). The filled black circles 
correspond to the positions of the atomic sites. 

are the creation and annihilation operators, respectively, 
of an electron at the site i with spin a. U is the strength 
of on-site Hubbard interaction. 

Decoupling of the interacting Hamiltonian: In or- 
der to determine the energy eigenvalues of the inter- 
acting model quantum system described by the tight- 
binding Hainiltonian given in Eq. [Jl first we decouple 
the interacting Hamiltonian using generalized Hartree- 
Fock approach, the so-called mean field approximation. 
In this approach, the full Hamiltonian is completely de- 
coupled into two parts. One is associated with the up- 
spin electrons, while the other is related to the down- 
spin electrons with their modified site energies. For 
up and down spin Hamiltonians, the modified site en- 
ergies are expressed in the form, eL = ei-\ -f U{ni_i) and 

^4 ^ ^il + U{ni^), where ni^ = cj^Cia is the number 
operator. With these site energies, the full Hamiltonian 
(Eq. [ij can be written in the decoupled form as. 
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where, Hf and H^ correspond to the effective tight- 
binding Hamiltonians for the up and down spin electrons, 
respectively. The last term is a constant term which pro- 
vides an energy shift in the total energy. 



Self consistent procedure: With these decoupled 
Hamiltonians (H-f and H^) of up and down spin elec- 
trons, we start our self consistent procedure considering 
initial guess values of (rii^) and {riii}. For these initial set 
of values of (rii-f) and (^14), we numerically diagonalize 
the up and down spin Hamiltonians. Then we calculate 
a new set of values of (riif) and (ri^). These steps arc 
repeated until a self consistent solution is achieved. 

Calculation of ground state energy: Using the self 
consistent solution, the ground state energy Eq for a par- 
ticular filling at absolute zero temperature (T = K) can 
be determined by taking the sum of individual states up 
to Fermi energy (Ep) for both up and down spins. Thus, 
we can write the final form of ground state energy as, 
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(3) 



where, the index n runs for the states upto the Fermi 
level. Enf (Eni) is the single particle energy eigenvalue 
for n-th eigenstate obtained by diagonalizing the Hamil- 
tonian H^ (Hi). 

Calculation of persistent current: At absolute zero 
temperature, total persistent current of the system is ob- 
tained from the expression. 
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where, Eq{(J)) is the ground state energy. 

Calculation of Drude weight: The Drude weight for 
the ring can be calculated through the relation, 

N f d^Eojcj)) 

where, N gives total number of atomic sites in the ring. 
Kohn^ has shown that for an insulating system D de- 
cays exponentially to zero, while it becomes finite for a 
conducting system. 

Determination of low-field magnetic susceptibil- 
ity: The general expression of magnetic susceptibility x 
at any flux tp is written in the form. 
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Evaluating the sign of x(0) "^c can able to predict 
whether the current is paramagnetic or diamagnetic in 
nature. Here we will determine x(0) only in the limit 
(/)—>■ 0, since we are interested to know the magnetic 
response in the low- field limit. 



B. Numerical results and discussion 

In this sub-section throughout our numerical work we 
set the nearest-neighbor hopping strength i = — 1 and 
second-neighbor hopping strength ti = —0.7. Energy 
scale is measured in unit of t. 



1. Perfect Hubbard Rings Described with NNH Integral 

For perfect rings we choose e^ = e^^ = for all i and 
since here we consider the rings with only NNH integral, 
the second-neighbor hopping strength ti is fixed at zero. 

Energy-flux characteristics: To explain the relevant 
features of magnetic response we begin with the energy- 
fiux characteristics. As illustrative examples, in Fig. [5] 
we plot the ground state energy levels as a function of 
magnetic flux (/) for some typical mesoscopic rings in the 
half-filled case, where (a) and (b) correspond to iV = 5 
and 6, respectively. The red curves represent the en- 
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FIG. 2: (Color online). Ground state energy levels as a func- 
tion of flux (j) for some typical mesoscopic rings in half-filled 
case. The red, green and blue curves correspond to [/ = 0, 1 
and 2, respectively, (a) A'' = 5 and (b) N — 6. 

crgy levels for the non-interacting ([/ = 0) rings, while 
the green and blue lines correspond to the energy levels 
for the interacting rings where the electronic correlation 
strength U is fixed to 1 and 2, respectively. From the 
spectra it is observed that the ground state energy level 
shifts towards the positive energy and it becomes much 
fiatter with the increase of the correlation strength U. 
Both for the two different ring sizes {N = 5 and 6) the 
ground state energy levels vary periodically with AB flux 
(j), but a significant difference is observed in their period- 
icities depending on the oddncss and evenness of the ring 
size N. For A^ = 6 (even), the energy levels show conven- 
tional (/)o (= 1, in our chosen unit system c = e = h = 1) 
fiux-quantum periodicity. On the other hand, the period 
becomes half i.e., 0o/2 for A^ = 5 (odd). This (/)o/2 peri- 
odicity disappears as long as the filling is considered away 
from the half-filling. At the same time, it also vanishes if 
impurities are introduced in the system, even if the ring 
is half-filled with odd A^. Therefore, 0o/2 periodicity is 




a special feature for odd half-filled perfect rings irrespec- 
tive of the Hubbard strength [/, while for all other cases 
traditional (j)o periodicity is obtained. 

To judge the accuracy of the mean-field calculations 
in our ring geometry, in Fig. [3] we show the variation of 
lowest energy levels where the eigenenergies are deter- 
mined through exact diagonalization of the full many- 
body Hamiltonian for the identical rings as given in 
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FIG. 3: (Color online). Ground state energy levels as a func- 
tion of flux (j) for some typical mesoscopic rings in half-filled 
case, where eigenenergies are determined through exact diago- 
nalization of the full many-body Hamiltonian. The red, green 
and blue curves correspond to {/ = 0, 1 and 2, respectively. 
(a) iV = 5 and (b) A^ = 6. 

Fig. [21 considering the same parameter values. Com- 
paring the results presented in Figs. [2] and [3l we see that 
the mean-field results agree very well with the exact di- 
agonalization method. Thus we can safely use mean-field 
approach to study magnetic response in our geometry. 

Current-flux characteristics: Following the above 
energy-flux characteristics now we describe the behavior 
of persistent current in mesoscopic Hubbard rings. As 
representative examples, in Fig. 2] we display the varia- 
tion of persistent currents as a function of flux (/> for some 
typical single-channel mesoscopic rings in the half-filled 
case, where (a) and (b) correspond to A'' = 15 and 20, 
respectively. The red, green and blue curves in Fig. Hl^a) 
correspond to the currents for [/ = 0, 1.5 and 2, re- 
spectively, while these curves in Fig. Hl^b) represent the 
currents for [/ = 0, 1 and 1.5, respectively. In the ab- 
sence of any e-e interaction ([/ = 0), persistent current 
shows saw-tooth like nature as a function of flux (f> with 
sharp transitions at n(j>o/2 (red line of Fig. IHa)) or ji(J>q 
(red line of Fig. Sl^b)), where n being an integer, de- 
pending on whether N is odd or even. The saw-tooth 



like behavior disappears as long as the electronic cor- 
relation is introduced into the system. This is clearly 
observed from the green and blue curves of Fig. SI Ad- 
ditionally, in the presence of U, the current amplitude 
gets suppressed compared to the current amplitude in 
the non-interacting case, and it decreases gradually with 
increasing U. This provides the lowering of electron mo- 
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FIG. 4: (Color online). Persistent current as a function of 
fiux for single-channel mesoscopic rings in half-filled case. 
(a) A'' = 15. The red, green and blue curves correspond to 
U = 0, 1.5 and 2, respectively, (b) A*' = 20. The red, green 
and blue curves correspond to [/ = 0, 1 and 1.5, respectively. 

bility with the rise of U and the reason behind this can 
be much better understood from our forthcoming discus- 
sion. Both for two different rings with sizes A^ = 15 (odd) 
and 20 (even), persistent currents vary periodically with 
AB flux (p showing different periodicities, following the 
energy-flux characteristics. For iV = 15, current shows 
(^o/2 flux-quantum periodicity, while for the other case 
(N = 20), current exhibits 0o flux-quantum periodicity. 

Variation of electronic mobility-Drude vireight: To 

reveal the conducting properties of Hubbard rings, we 
study the variation of Drude weight D for these systems. 
Drude weight can be calculated by using Eq. [5j Finite 
value of D predicts the metallic phase, while for the in- 
sulating phase it drops exponentially to zero^. 

As illustrative examples, in Fig. [S] we show the varia- 
tion of Drude weight Z? as a function of electronic correla- 
tion strength U for some typical single-channel Hubbard 
rings. In Fig.ISJa) the results are shown for three different 
half-filled rings, where the red, green and blue lines corre- 
spond to the rings with N = 10, 30 and 50, respectively. 
From the curves it is evident that for smaller values of t/, 
the half-filled rings show finite value of D which reveals 
that they are in the metallic phase. On the other hand. 



D drops sharply to zero when U becomes high. Thus 
the rings become insulating when U is quite large. The 
results for the non-half filled case are shown in Fig. EJb), 
where we fix the ring size A'^ = 20 and vary the electron 
filling. The red, green and blue curves represent iVg = 10, 
14 and 18, respectively, where Ne gives the total number 
of electrons in the ring. For these three choices of Ne, 
the ring is always less than half-filled (since N^ < N) 
and the ring is in the conducting phase irrespective of 



half-filled rings, ASD is less than unity for small value of 
U, while it reaches to unity when U is large. This be- 
havior is clearly shown by the curves given in Fig. [HJa), 
where the red, green and blue lines correspond to ASDs 
for the half-filled rings with TV = 10, 30 and 50, respec- 
tively. Thus, for low U there is some finite probability of 
getting two opposite spin electrons in a same site which 
allows electrons to move along the ring and the metallic 
phase is obtained. But for large U, ASD reaches to unity 
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FIG. 5: (Color online). Drude weight as a function of Hub- 
bard interaction strength U for single-channel niesoscopic 
rings, (a) Half-filled case. The red, green and blue curves 
correspond to N = 10, 30 and 50, respectively, (b) Non- 
half-filled case with A'^ = 20. The red, green and blue curves 
correspond to A^e = 10, 14 and 18, respectively. 



FIG. 6: (Color online). Average spin density (ASD) as a 
function of Hubbard interaction strength U for single-channel 
niesoscopic rings, (a) Half-filled case. The red, green and 
blue curves correspond to A = 10, 30 and 50, respectively, 
(b) Non-half-filled case with N = 20. The red, green and 
blue curves correspond to Ne = 10, 14 and 18, respectively. 



the correlation strength U. Now we try to justify the 
dependence of the Hubbard strength U on the electronic 
mobility for these different fillings. To understand the 
effect of U on electron mobility here we measure a quan- 
tity called 'average spin density' (ASD) which is defined 
by the factor ^^ \{ni-f — ni^)\/N. The integer i is the site 
index and it runs from 1 to N . By calculating ASD we 
can estimate the occupation probability of electrons in 
the ring and it supports us to explain whether the ring 
lies in the metallic phase or in the insulating one. For 
the rings those are below half-filled, ASD is always less 
than unity irrespective of the value of U as shown by 
the curves in Fig. [BJb). It reveals that for these systems, 
ground state always supports an empty site and elec- 
tron can move along the ring avoiding double occupancy 
of two different spin electrons at any site i in the pres- 
ence of e-e correlation which provides the metallic phase 
{D > 0). For a fixed ring size and a particular strength 
of t/, the ASD increases as the filling is increased towards 
half-filling which is noticed by comparing the three dif- 
ferent curves in Fig. IHl^b). On the other hand, in the 



which means that each site is singly occupied either by 
an up or down spin electron with probability 1. In this 
case ground state docs not support any empty site and 
due to strong repulsive e-e correlation one electron sit- 
ting in a site does not allow to come other electron with 
opposite spin from the neighboring site which provides 
the insulating phase [D = 0). The situation is some- 
what analogous to Mott localization in one-dimensional 
infinite lattices. In perfect Hubbard rings the conducting 
nature has been studied exactly quite a long ago using 
the ansatz of Bethe by Shastry and Sutherland^. They 
have calculated charge stiffness constant {Dc) and have 
predicted that Dc goes to zero as the system approaches 
towards half-filling for any non-zero value of U . Our nu- 
merical results clearly justify their findings. 

Low-field magnetic susceptibility: Now, we discuss 
the variation of low-field magnetic susceptibility which 
can be calculated from Eq. [Blby setting — >■ 0. With the 
help of this parameter we can justify whether the cur- 
rent is paramagnetic (-|-ve slope) or diamagnetic (— ve 



slope) in nature. For our illustrative purposes, in Fig. [7] 
we show the variation of low-field magnetic susceptibil- 
ity with system size N for some typical single-channel 
mesoscopic rings in the half-filled case. Figure ^a.) cor- 
respond to the variation of low-field magnetic suscepti- 
bility for the non-interacting {U = 0) rings, where the 
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FIG. 7: (Color online). Low-field magnetic susceptibility as a 
function of system size A'^ for single-channel mesoscopic rings 
in half-filled case, (a) U = 0. A*' is an odd {2n+ 1) or an even 
(2?i + 2) number, where n is an integer, (b) U — 1. A'' is an 
even number obeying the relation N — 4n + 2. (c) U = 1. N 
is an even number satisfying the relation A'^ = 4n. (d) U — 1. 
N is an odd number following the relation A^ = 2n -|- 1. 

ring size can by anything i.e., either odd, following the 
relation iV = 2n + 1 (n is an integer), or even, obey- 
ing the expression N = 2n + 2. ft is observed that 
both for odd and even TV, low-field current exhibits dia- 
magnetic nature. The behavior of the low-field currents 
changes significantly when the e-e interaction is taken 
into account. Depending on the ring size N, the sign 
becomes -l-ve and — ve as shown by the curves given in 
Figs. [7{b)-(d). For the interacting rings where the rela- 
tion A^ = 4n-|-2 is satisfied, the low-field current becomes 
diamagnetic (Fig.[71Jb)). The sign becomes paramagnetic 
when N = An (Fig. Ulc)) a.nd N = 2n + 1 (Fig. [ltd)). 
Thus, in brief, we say that for non-interacting half-filled 
rings low-field current exhibits diamagnetic response ir- 
respective of N i.e., whether N is odd or even. For the 
interacting half-filled rings with odd iV, low-field current 
provides only the paramagnetic behavior, while for even 
A'^, depending on the particular value of N, the response 
becomes either diamagnetic or paramagnetic. These na- 
tures of low-field currents change for the cases of other 
electron fillings. Hence, it can be emphasized that the be- 
havior of the low-field currents is highly sensitive on the 
Hubbard correlation, electron filling, evenness and odd- 



ness of N, etc. The behavior of zero-field magnetic sus- 
ceptibility in Hubbard rings has been studied extensively 
quite a long back using the Bethe ansatz by Shibap^. 
In this work, he has studied magnetic susceptibility per 
electron as functions of electron filling and Hubbard cor- 
relation strength and provided several interesting results. 
From his findings we can clearly justify our presented re- 
sults. 



2. Disordered Hubbard Rings Described with NNH and 
SNH Integrals 

Now, we explore the combined effect of electron- 
electron correlation and second-neighbor hopping (SNH) 
integral on persistent current in disordered mesoscopic 
rings. 

To get a disordered ring, we choose site energies (eif 
and Cii) randomly from a "Box" distribution function of 
width W. As the site energies are chosen randomly it 
is needed to consider the average over a large number of 
disordered configurations (from the stand point of statis- 
tical average). Here, we determine the currents by taking 
the average over 50 random disordered configuration in 
each case to achieve much accurate results. 

As illustrative examples, in Fig. [5] we display the varia- 
tion of persistent currents for some single-channel meso- 
scopic rings considering 1/3 electron filling. In (a) the 
results are given for the rings characterized by the NNH 
integral model. The red curve represents the current for 
the ordered {W = 0) non-interacting {U = 0) ring. It 
shows saw-tooth like nature with AB flux providing (Jjq 
flux-quantum periodicity. The situation becomes com- 
pletely different when impurities are introduced in the 
ring as clearly seen by the other two colored curves. The 
green curve represents the current for the case only when 
impurities are considered but the effect of Hubbard inter- 
action is not taken into account. It varies continuously 
with (f> and gets much reduced amplitude, even an or- 
der of magnitude, compared to the perfect case. This is 
due to the localization of the energy eigenstates in the 
presence of impurity, which is the so-called Anderson lo- 
calization. Hence, a large difference exists between the 
current amplitudes of an ordered and disordered non- 
interacting rings and it was the main controversial is- 
sue among the theoretical and experimental predictions. 
Experimental results suggest that the measured current 
amplitude is quite comparable to the theoretically esti- 
mated current amplitude in a perfect system. To remove 
this controversy, as a first attempt, we include the effect 
of Hubbard interaction in the disordered ring described 
by the NNH model. The result is shown by the blue 
curve where U is fixed at 0.5. It is observed that the 
current amplitude gets increased compared to the non- 
interacting disordered ring, though the increment is too 
small. Not only that the enhancement can take place 
only for small values of U, while for large enough U the 
current amplitude rather decreases. This phenomenon 



can be explained as follows. For the non-interacting dis- 
ordered ring the probability of getting two opposite spin 
electrons becomes higher at the atomic sites where the 
site energies are lower than the other sites since the elec- 
trons get pinned at the lower site energies to minimize 
the ground state energy, and this pinning of electrons 
becomes increased with the rise of impurity strength W. 
As a result the mobility of electrons and hence the cur- 
rent amplitude gets reduced with the increase of impurity 
strength W. Now, if we introduce electronic correlation 
in the system then it tries to depin two opposite spin 
electrons those are situated together due to the Coulomb 
repulsion. Therefore, the electronic mobility is enhanced 
which provides quite larger current amplitude. But, for 




FIG. 8: (Color online). Persistent current as a function of flux 
(j> for single-channel mesoscopic rings with Ai" = 15 considering 
1/3 electron flUing. (a) Rings with only NNH integral. The 
red line corresponds to the ordered non-interacting ring, while 
the green and blue lines correspond to the disordered {W = 2) 
rings with (7 = and 0.5, respectively, (b) Rings with NNH 
and SNH integrals. The red line represents the ordered non- 
interacting ring, whereas the green and blue line correspond to 
the disordered {W = 2) rings with U = and 1.5, respectively. 



is clearly observed from Fig. [IJb). The red curve refers 
to the current for the perfect {W = 0) non-interacting 
([/ = 0) ring and it achieves much higher amplitude com- 
pared to the NNH model (see red curve of Fig. [5{a)). 
This additional contribution comes from the SNH inte- 
gral since it allows electrons to hop further. The main 
focus of this sub-section is to interpret the combined ef- 
fect of SNH integral and Hubbard correlation on the en- 
hancement of persistent current in disordered ring. To 
do this first we narrate the effect of SNH integral in dis- 
ordered non-interacting ring. The nature of the current 
for this particular case is shown by the green curve of 
Fig. EKb). It shows that the current amplitude gets re- 
duced compared to the perfect case (red line), which is 
expected, but the reduction of the current amplitude is 
very small than the NNH integral model (see green curve 
of Fig.mja)). This is due the fact that the SNH integral 
tries to delocalize the electronic states, and therefore, 
the mobility of the electrons is enriched. The situation 
becomes more interesting when we include the effect of 
Hubbard interaction. The behavior of the current in the 
presence of interaction is plotted by the blue curve of 
Fig.[5Jb) where we fix [/ = 1.5. Very interestingly we see 
that the current amplitude is enhanced moderately and 
quite comparable to that of the perfect ring. Therefore, 
it can be predicted that the presence of SNH integral 
and Hubbard interaction can provide a persistent cur- 
rent which inay be comparable to the measured current 
amplitudes. In the above analysis we consider the effect 
of only SNH integral in addition to the NNH model, and, 
illustrate how such a higher order hopping integral leads 
to an important role on the enhancement of current am- 
plitude in presence of Hubbard correlation for disordered 
rings. Instead of considering only the SNH integral we 
can also take the contributions from all possible higher 
order hopping integrals with reduced hopping strengths. 
Since the strengths of other higher order hopping inte- 
grals are too small, the contributions from these factors 
are reasonably small and they will not provide any sig- 
nificant change in the current amplitude. 



III. A HUBBARD CYLINDER IN ABSENCE OF 
SO INTERACTIONS 



large enough interaction strength, mobility of electrons 
gradually decreases due to the strong repulsive interac- 
tion. Accordingly, the current amplitude gradually de- 
creases with U. So, in short, we can say that within the 
nearest-neighbor hopping (NNH) model electron-electron 
interaction does not provide any significant contribution 
to enhance the current amplitude, and hence the contro- 
versy regarding the current amplitude still persists. 

To overcome this controversy, finally we make an at- 
tempt by incorporating the effect of second-neighbor hop- 
ping (SNH) integral in addition to the nearest-neighbor 
hopping (NNH) integral. With this modification a sig- 
nificant change in current amplitude takes place which 



In this section we extend our discussion for an inter- 
acting mesoscopic ring with finite width threaded by an 
AB fiux (p. Here also we ignore the effect of SO inter- 
action on magneto-transport properties like the previous 
section. 



A. Model and the Hamiltonian 

Let us start by referring to Fig.|9l where a small metal- 
lic cylinder is threaded by a magnetic fiux 0. The filled 
black circles correspond to the positions of the atomic 
sites in the cylinder. To predict the size of a cylinder we 



use two parameters N and M, where the 1st one (N) rep- 
resents total number of atomic sites in each circular ring 
and the other one (M) gives total number of identical 
circular rings. For the description of our model quantum 
system we use a tight-binding framework and in order to 
incorporate the effect of higher order hopping integrals to 
the Hamiltonian here we consider second-neighbor hop- 
ping (SNH) (shown by the red dashed line in Fig. ^ in 
addition to the nearest-neighbor hopping (NNH) of elec- 
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FIG. 9: (Color online). Schematic view of a ID mesoscopic 
cylinder penetrated by a magnetic flux (p. The red dashed 
line corresponds to the second-neighbor hopping integral and 
the flUed black circles represent the positions of the atomic 
sites. A persistent current / is established in the cylinder. 

trons. Considering both NNH and SNH integrals the TB 
Hamiltonian for the cylindrical system in Wannier basis 
looks in the form, 



■ij.c 



1,3,17 



h.c] + ^ id le*^''cj ,,■ ,^Ci+i,j+i,CT + h.c. 



(7) 



where, (i,j) represent the co-ordinate of a lattice site. 
The index i runs from 1 to M, while the integer j goes 
from 1 to N. eij,cr is the on-site energy of an electron 
at the site {i,j) of spin a (t,i). ti and td arc the NNH 
and SNH integrals, respectively. Due to the presence of 
magnetic flux (j) (measured in unit of the elementary flux 
quantum 0o = ch/e) , a phase factor 9i = 2tt(J)/N appears 
in the Hamiltonian when an electron hops longitudinally 
from one site to its neighboring site, and accordingly, a 
negative sign comes when the electron hops in the re- 
verse direction. 9d is the associated phase factor for the 
diagonal motion of an electron between two neighboring 
concentric rings. No phase factor appears when an elec- 
tron moves along the vertical direction which is set by 
proper choice of the gauge for the vector potential A as- 
sociated with the magnetic field B, and this choice makes 
the phase factors (6*/, 9d) identical to each other for the 
longitudinal and diagonal motions. Since the magnetic 
field corresponding to the AB flux (p does not penetrate 
anywhere of the surface of the cylinder, we ignore Zeeman 
term in the above tight-binding Hamiltonian (Eq. [7|). 



*J,cr 



and 



'ij,^ 



are the creation and annihilation oper- 



ators, respectively, of an electron at the site {i,j) with 
spin a. U is the on-site Hubbard interaction term. 



B. Theoretical formulation 

To calculate energy eigenvalues, persistent current and 
related issues here we follow exactly the same prescrip- 
tion which we illustrate in the earlier section (Sec. II) 
i.e.; Hartree-Fock mean- field approach. 



C. Numerical results and discussion 

Throughout the numerical analysis, in this sub-section, 
we set the nearest-neighbor hopping strength </ = — 1 and 
fix M = 2 i.e., cylinders with two identical rings. Energy 
scale is measured in the unit of ti . We describe the re- 
sults in three different parts. In the first part, we consider 
perfect cylinders with only nearest-neighbor hopping in- 
tegral. In the second part, disordered cylinders described 
with only NNH integral are considered. Finally, in the 
third part we discuss the effect of second-neighbor hop- 
ping (SNH) integral on the enhancement of persistent 
current in disordered cylinders. 



1. Perfect cylinders with NNH integral 

For perfect cylinders we choose e^ j.-j- = £i.j,i = for all 
(i, j). Since here we consider the cylinders described with 
NNH integral only, the second-neighbor hopping strength 
td is fixed to zero. 

Energy-flux characteristics: As illustrative examples, 
in Fig. [TO]wc show the variation of ground state energy 
levels as a function of magnetic flux (jj for some typical 
mesoscopic cylinders where N is fixed at 5 (odd N). In 
(a) the results are given for the quarterly-filled (N^ ~ 5) 
cylinders, while in (b) the curves correspond to the re- 
sults for the half-filled {N^ = 10) cylinders. The red, 
green and blue lines represent the ground state energy 
levels for U — 0, 0.5 and 1, respectively. It is observed 
that the ground state energy shows oscillatory behavior 
as a function of (/> and the energy increases as the elec- 
tronic correlation strength U gets increased. Most signif- 
icantly we see that the ground state energy levels provide 
two different types of periodicities depending on the elec- 
tron filling. At quarter-filling, ground state energy level 
gives (j)Q (= 1, since c = e = h=lm our chosen unit 
system) flux-quantum periodicity. On the other hand, at 
half-fllling it shows (/)o/2 flux-quantum periodicity. The 
situation becomes quite different when the total number 
of atomic sites N in individual rings is even. For our 
illustrative purposes in Fig. [TT] we plot the lowest en- 
ergy levels as a function of (j) for some typical mesoscopic 
cylinders considering iV = 8 (even TV). The curves of 
different colors correspond to the identical meaning as in 
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Fig. [TOl From the spectra given in Figs. [TlT a') (quarter- 
filled case) and (b) (half-filled case) it is clearly observed 
that the ground state energy levels vary periodically with 
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FIG. 10: (Color online). Ground state energy levels as a 
function of flux (j) for some perfect cylinders with N — 5 and 
M — 2. The red, green and blue curves correspond to f/ = 0, 
0.5 and 1, respectively, (a) Quarter-filled case and (b) Half- 
filled case. 
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FIG. 11: (Color online). Ground state energy levels as a 
function of flux (p for some perfect cylinders considering N — 8 
and Ad = 2. The red, green and blue curves correspond to 
[7 = 0, 0.5 and 1, respectively, (a) Quarter-filled case and (b) 
Half-filled case. 

AB flux (f) exhibiting only (j)o flux-quantum periodicity. 
Thus it can be emphasized that the appearance of half 
flux-quantum periodicity strongly depends on the elec- 
tron filling as well as on the oddness and evenness of the 
total number of atomic sites N in individual rings. Only 
for the half-filled cylinders with odd N, the lowest energy 



level gets (j)o/2 periodicity with flux (p. Now it is impor- 
tant to note that this half flux-quantum periodicity does 
not depend on the width (M) of the cylinder and also 
it is independent of the Hubbard correlation strength U. 
Hence, depending on the system size and filling of elec- 
trons variable periodicities are observed in the variation 
of lowest energy level. It may provide an important sig- 
nature in studying magnetic response in nano-scale loop 
geometries. 

Current-flux characteristics: In Fig. [T^] we display 
the current-flux characteristics for some impurity free 
mesoscopic cylinders considering M = 2. In (a) the 
results are given for the half-filled case where we set 
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FIG. 12: (Color online). Persistent current as a function of 
fiux (j) for some ordered mesoscopic cylinders considering M = 
2. (a) Half-filled case with A'^ — 15. The red, green and blue 
curves correspond to U — 0, 1.5 and 2, respectively, (b) 
Quarter-filled case with A*' = 20. The red, green and blue 
curves correspond to f/ = 0, 2 and 3, respectively. 

N — 15. The red line corresponds to the current for 
the non-interacting [U = 0) case, while the green and 
blue lines represent the currents when U — 1.5 and 2, 
respectively. From the curves we notice that the cur- 
rent amplitude gradually decreases with the increase of 
electronic correlation strength U. The reason is that at 
half-filling each site is occupied by at least one electron 
of up spin or down spin, and the placing of a second 
electron of opposite spin needs more energy due to the 
repulsive effect of U. Thus conduction becomes difficult 
as it requires more energy when an electron hops from 
its own site and situates at the neighboring site. Now 
both for the non-interacting and interacting cases, cur- 
rent shows half flux-quantum periodicity as a function a 
(j) obeying the energy-flux characteristics since here we 
choose odd N {N = 15). The behavior of the persis- 
tent currents for even N is shown in (b) where we set 
A'^ = 20. The currents are drawn for the quarter-filled 
case i.e., N,, = 20, where the red, green and blue curves 
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correspond to C/ = 0, 2 and 3, respectively. The reduc- 
tion of current amplitude with the increase of Hubbard 
interaction strength is also observed for this quarter-filled 
case, similar to the case of half-filled as described earlier. 
But the point is that at quarter-filling, the reduction of 
current amplitude is much smaller compared to the half- 
filled situation. This is quite obvious in the sense that 
at less than half-filling 'empty' lattice sites are available 
where electrons can hop easily without any cost of ex- 
tra energy and the conduction becomes much easier than 
the half-filled situation. In this quarter-filled case, per- 
sistent currents provide only (J)q flux-quantum periodicity 
following the E-(j) diagram. From these current-flux char- 
acteristics it can be concluded that for 'ordered ' cylinders 
current amplitude always decreases with the enhancement 
in Hubbard correlation strength U . 



2. Disordered cylinders with NNH integral 

In order to describe the effect of impurities on elec- 
tron transport now we focus our attention on the re- 
sults of some typical disordered cylinders described with 
NNH integral. Here we consider the diagonal disordered 
cylinders i.e., impurities are introduced only at the site 

-15.70 I 



^ -15.85 



16.00 




-1 -0.5 O 0.5 1 



-28.70 



^ -28.95 




-29.20 



-1 -0.5 O 0.5 1 



FIG. 13: (Color online). Ground state energy level as a func- 
tion of flux (j) for half-flUed disordered mesoscopic cylinders 
(M = 2) considering U = 1 and W = 2. (a) N = 5 and (b) 

iV = 8. 

energies without disturbing the hopping integrals. The 
site energies in each concentric ring are chosen from a 
correlated distribution function which looks in the form, 



ej,t = ^j-l = WcosijXn) 



(8) 



where, W is the impurity strength. A is an irrational 
number and we choose A — (l + \/5)/2, for the sake of our 
illustration. Setting A = 0, we get back the pure system 
with uniform site energy W. Now, instead of considering 



site energies from a correlated distribution function, as 
mentioned above in Eq. [51 we can also take them ran- 
domly from a "Box" distribution function of width W. 
But in the later case we have to take the average over 
a large number of disordered configurations (from the 
stand point of statistical average) and since it is really a 
difficult task in the aspect of numerical computation we 
select the other option. Not only that in the averaging 
process several mesoscopic phenomena may disappear. 
Therefore, the averaging process is an important issue in 
low-dimensional systems. 

In presence of disorder, energy levels get modified sig- 
nificantly. For our illustrative purposes in Fig.[T3]we plot 
ground state energy levels as a function of magnetic fiux 
(j) for some disordered mesoscopic cylinders when they 
are half-filled. The Hubbard interaction strength U is 
set at 1 and the impurity strength W is fixed to 2. In (a) 
the ground state energy level is shown for a cylinder with 
N = 5 (odd), while in (b) it is presented for a cylinder 
taking A^ = 8 (even). Quite interestingly we see that for 
the cylinder with odd iV, the half flux-quantum period- 
icity of the lowest energy level disappears in the presence 
of impurity and it provides conventional 0o periodicity. 
Hence, for cylinders with odd N, 0o/2 flux-quantum pe- 
riodicity will be observed only when they are free from 
any impurity. For the disordered cylinder with even N 
{N = 8), the lowest energy level as usual provides 00 pe- 
riodicity similar to the impurity free cylinders containing 
even N. Apart from this periodic nature, impurities play 
another significant role in the determination of the slope 
of the energy levels. The slope of the lowest energy level 
decreases significantly compared to the perfect case, and 
therefore, a prominent change in current amplitude also 
takes place. 

To justify the above facts, in Fig. [TJ] we present the 
variations of persistent currents with AB flux for a half- 
filled mesoscopic cylinder, described in the framework of 
NNH model, considering iV = 15 and M = 2. The red 
curve represents the current for the ordered {W = 0) 
non-interacting (U ~ 0) cylinder. It shows saw-tooth like 
nature with flux providing 4>o/2 flux-quantum period- 
icity. The situation becomes completely different when 
impurities are introduced in the cylinder as seen by the 
other two curves. The green curve represents the current 
for the case only when impurities are considered but the 
effect of electronic correlation is not taken into account. 
It shows a continuous like nature with (f>Q flux-quantum 
periodicity. The most important observation is that the 
current amplitude gets reduced enormously, even an or- 
der of magnitude, compared to the perfect cylinder. This 
is due to the localization of the energy eigenstates in the 
presence of impurity, which is the so-called Anderson lo- 
calization. Hence, a large difference exists in the current 
amplitudes of an ordered and disordered non-interacting 
cylinders and it was the main controversial issue among 
the theoretical and experimental predictions. Experi- 
mental verifications suggest that the measured current 
amplitude is quite comparable to the theoretical current 
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amplitude obtained in a perfect system. To remove this 
controversy, as a first attempt, we include the effect of e- 
e correlation in the disordered cylinder described by the 
NNH model. The result is shown by the blue curve where 
U is fixed at 1.5. It is observed that the current am- 
plitude gets increased compared to the non-interacting 
disordered cylinder, though the increment is too small. 
Not only that the enhancement can take place only for 
small values of U, while for large enough U the current 
amplitude rather decreases. This phenomenon can be im- 
plemented as follows. For the non-interacting disordered 
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FIG. 14: (Color online). Persistent current as a function of 
flux for a half- filled mesoscopic cylinder considering TV = 15 
and M = 2. The red line corresponds to the ordered case 
when U = 0, whereas the green and blue lines correspond to 
the disordered case (W = 2) when U = and 1.5, respectively. 



cylinder the probability of getting two opposite spin elec- 
trons becomes higher at the atomic sites where the site 
energies are lower than the other sites since the electrons 
get pinned at the lower site energies to minimize the 
ground state energy, and this pinning of electrons be- 
comes increased with the rise of impurity strength W. 
As a result the mobility of electrons and hence the cur- 
rent amplitude gets reduced with the increase of impurity 
strength W. Now, if we introduce electronic correlation 
in the system then it tries to depin two opposite spin 
electrons those are situated together due to the Coulomb 
repulsion. Therefore, the electronic mobility is enhanced 
which provides larger current amplitude. But, for large 
enough interaction strength, no electron can able to hop 
from one site to other at the half-filling since then each 
site is occupied either by an up or down spin electron 
which does not allow other electron of opposite spin due 
to the repulsive term U. Accordingly, the current ampli- 
tude gradually decreases with U. On the other hand, at 
less than half-filling though there is some finite probabil- 
ity to hop an electron from one site to the other available 
'empty' site but still it is very small. So, in brief, we 
can say that within the nearest-neighbor hopping (NNH) 
approximation electron-electron interaction does not pro- 
vide any significant contribution to enhance the current 
amplitude, and hence the controversy regarding the cur- 
rent amplitude still persists. 



3. Disordered cylinders with NNH and SNH integrals 

To overcome the existing situation regarding the cur- 
rent amplitude, in this sub-section, finally we make an at- 
tempt by incorporating the effect of second-neighbor hop- 
ping (SNH) integral in addition to the nearest-neighbor 
hopping (NNH) integral. 

A significant change in current amplitude takes place 
when we include the contribution of second-neighbor hop- 
ping (SNH) integral in addition to the NNH integral. As 
representative examples, in Fig. [15] we plot the current- 
flux characteristics for a half-filled mesoscopic cylinder 
considering iV = 15 and M = 2. The black, magenta 
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FIG. 15: (Color online). Persistent current as a function of 
flux <f) for a half- filled mesoscopic cylinder taking N = 15 and 
M = 2 in the presence of NNH and SNH integrals. The black 
line corresponds to the ordered case when U = 0, whereas 
the magenta and gold lines correspond to the disordered case 
{W = 2) when U — and 1.5, respectively. Here SNH integral 
is fixed at —0.6. The currents shown by the red, green and 
blue lines for the ring described with NNH model (identical 
to Fig. [14} are re-plotted to judge the effect of SNH integral 
over NNH model much clearly. 

and gold lines correspond to the results in the presence 
of SNH integral, while the other three colored curves 
(red, green and blue) represent the currents in the ab- 
sence of SNH integral. Here we choose td = —0.6. The 
black curve refers to the persistent current for the perfect 
{W = 0) non-interacting {U — 0) cylinder and it achieves 
much higher amplitude compared to the NNH model (red 
curve). This additional contribution comes from the SNH 
integral since it allows electrons to hop further. In addi- 
tion it is also noticed that the current varies periodically 
with (j> providing 0o flux-quantum periodicity, instead of 
0o/2 as in the case of NNH integral model (red curve). 
Thus, it can be emphasized that (/)o/2 periodicity will be 
observed only when the cylinder is (a) free from impurity, 
(b) half-filled, (c) made with odd N, and (d) described 
by the nearest-neighbor hopping model. The main focus 
of this sub-section is to interpret the combined effect of 
SNH integral and electron-electron correlation on the en- 
hancement of persistent current amplitude in disordered 
cylinder. To do this first we narrate the effect of SNH 
integral in disordered non-interacting cylinder. The na- 
ture of the current for this particular case is shown by the 
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magenta curve of Fig. [15] It shows that the current am- 
pUtude gets reduced compared to the perfect case (black 
line) , which is expected, but the reduction of the current 
amphtude is quite small than the NNH integral model. 
This is due the fact that the SNH integral tries to de- 
localize the electronic states, and therefore, the mobility 
of the electrons is enriched. The situation becomes more 
interesting when we include the effect of Hubbard inter- 
action. The behavior of the current in the presence of 
interaction is plotted by the gold curve of Fig. [T5] where 
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FIG. 16: (Color online). Persistent current as a function of 
flirx (j) for a half-filled mesoscopic cylinder taking A'' = 15 and 
M = 2 in the presence of NNH and SNH integrals. The black 
line corresponds to the ordered case when U = 0, whereas 
the magenta and gold lines correspond to the disordered case 
{W = 2) when U — and 1.5, respectively. Here SNH integral 
is fixed at —0.8. The currents shown by the red, green and 
blue Unes for the ring described with NNH model (identical 
to Fig. I14p are re-plotted to judge the effect of SNH integral 
over NNH model much clearly. 

we fix [/ ~ 1.5. Very interestingly we see that the current 
amplitude is enhanced significantly and quite comparable 
to that of the perfect cylinder. 

For better clarity of the results discussed above, in 
Fig. [Tni we also present the similar feature of persistent 
current for other hopping strength of SNH integral. Here 
we set td = —0.8. From these curves we see that the 
current amplitude gets enhanced more as we increase the 
SNH strength. 

Thus, it can be emphasized that the presence of higher 
order hopping integrals and electron-electron correlation 
may provide a persistent current which can be compara- 
ble to the measured current amplitudes. Throughout the 
above analysis we set the width of the cylinders at a fixed 
value (M = 2), for the sake of our illustration. But, all 
these results are also valid for cylinders of larger widths. 



IV. A MESOSCOPIC RING WITH RASHBA 
AND DRESSELHAUS SO INTERACTIONS 

Finally, in this section we address the magneto- 
transport properties in a mesoscopic ring, threaded by 
an AB fiux </>, in the presence of Rashba and Dressclhaus 
SO interactions. We establish that the presence of SO 



interaction, in general, leads to an enhanced amplitude 
of the persistent current. This is another one approach 
through which we can justify the appearance of larger 
current amplitude in a disordered ring. For this discus- 
sion we neglect the effect of electron-electron interaction. 



Model, TB Hamiltonian and the theoretical 
formulation 



The schematic view of a mesoscopic ring subjected to 
an AB flux (jj (measured in unit of the elementary fiux 
quantum (J)q = ch/e) is shown in Fig. 1171 Within a TB 




FIG. 17: (Color online). A mesoscopic ring threaded by an 
AB flux (f). 

framework the Hamiltonian for such an iV-site ring is^Ti^ 
(and references therein), 



H — Hq + Hg 



Here, 



^0 = H (^UoCr. 






{clte"'cr,+i + h.c: 



and, 



where. 



Hso 






^n^ t .«! 






h.c. 



(9) 



(10) 



(11) 



= Urso {(Tx cos (/3„,„+l -I- (Ty sin (pn,n+l ) 

- itoso (cTy cos(pn,n+i + ^x sin(p„^„+i) . 



(12) 



n=l, 2, ...,iVis the site index along the azimuthal 
direction ip of the ring. The other factors in Egs. 1101 and 
[Tulare as follows. 



Here eo is the site energy of each atomic site of 
the ring, t is the nearest-neighbor hopping integral 
and Q = 27t<j>/N is the phase factor due to the AB 
flux (f> threaded by the ring, tjiso and toso are the 
isotropic nearest-neighbor transfer integrals which 
measure the strengths of Rashba and Dresselhaus SO 
couplings, respectively, and (pn,n+i = i^n + ^n+i) /'2, 
where ipn = 27r(n — ^)/N. ax and a-y are the Pauli 
spin matrices, cjjg. (cna) is the creation (annihilation) 
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operator of an electron at the site n with spin a (t,i)- 
Throughout the numerical analysis, in this sub-section, 
we choose t — 1 and measure the SO coupling strength 
in unit of t. 

At absolute zero temperature (T = OK), the persis- 
tent current in the ring described with fixed number of 
electrons Ne is determined by. 



/((/>) = -, 



' dcj) 



(13) 



where, Eq{(J}) is the ground state energy. Wc compute 
this quantity by exactly diagonalizing the TB Hamilto- 
nian (Eq. [5]) to understand unambiguously the role of the 
RSOI interaction alone on persistent current. 



B. Numerical results and discussion 



Energy-flux characteristics: Before presenting the 
results for /(^), to make the present communication a 
self contained study, we first take a look at the energy 
spectrum of both an ordered and a disordered ring with 




-0.5 



0.5 
>Flux 




FIG. 18: (Color online). E-(j) curves of a 12-site ring, where 
the 1st and 2nd columns correspond to the results for the 
ordered {W = 0) and disordered {W = 1) cases, respectively. 
The red, green and blue lines correspond to Irso = toso = 0; 
tflso = 1, toao = and Irso = 1, toao = 0.5, respectively. 

and without the SO interactions, as the fiux through the 
ring is varied. In Fig. [18] the flux dependent spectra are 
shown for a 12-site ordered ring and a randomly disor- 
dered one (with diagonal disorder) in the left and the 
right columns respectively. Clearly, disorder destroys the 
band crossings observed in the ordered case. The pres- 
ence of the RSOI and the DSOI also lifts the degeneracy 
and opens up gaps towards the edges of the spectrum. 



Enhancement of persistent current: 

• An ordered ring: In Fig. [TO] we examine the effect of the 
RSOI on the persistent current of an ordered ring with 80 
sites. The DSOI is set equal to zero. We have examined 
both the non-half-filled and half-filled band cases, but 
present results for the latter only to save space. With 
increasing strength of the RSOI the persistent current 




FIG. 19: (Color online). Current-flux characteristics of a 80- 
site ordered {W = 0) half-filled ring for different values of 
tRso when tDso is set at 0. 

exhibits a trend of an increase in its amplitude. Local 
phase reversals take place together with the appearance 
of kinks in the current-fiux diagrams which are however, 
not unexpected even without the RSOI, and are results of 
the band crossings observed in the spectra of such rings. 
The amplitude of the persistent current at a specific value 
of the magnetic flux is of course not predictable in any 
simple manner, and is found to be highly sensitive to 
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FIG. 20: (Color online). Persistent current at a particular AB 
flux (0 — 0.25) as a function Rashba SO interaction strength 
for an ordered {W = 0) half-filled ring with N — 60 when 
toso is set to zero. 

the number of electrons N^ (i.e., the fllling factor). Is- 
sues related to the dependence of the persistent current 
on the filling factor have been elaborately discussed by 
Splettstoesser et al^. 

The persistent current in an ordered ring also exhibits 
interesting oscillations in its amplitude as the RSOI is 
varied keeping the magnetic fiux fixed at a particular 
value. The oscillations persist irrespective of the band- 
filling factor Ne, with or without the presence of the 
DSOI. In Fig. [5n]the oscillating nature of persistent cur- 
rent is presented for a 60-site ordered ring in the half- 
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filled band case when </> is set at (/>o/4. The current ex- 
hibits oscillations with growing amplitude as the strength 
of the RSOI is increased. 

• A disordered ring: We now present the results for a 
disordered ring of 80 sites in Fig. [21] Disorder is intro- 
duced via a random distribution (width M^ = 2) of the 
values of the on-site potentials (diagonal disorder), and 
results averaged over sixty disorder configurations have 
been presented. The DSOI remains zero. Without any 
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FIG. 21: (Color online). Current-flux characteristics of a 80- 
site disordered {W = 2) half-fllled ring for different values of 
tRao when toso is fixed at 0. 

spin-orbit interaction, disorder completely suppresses the 
persistent current (an effect of the localization of the elec- 
tronic states in the ring), as it is observed in Fig. [2T] 
(black curve). With the introduction of the RSOI, the 
current starts increasing, and for tuso — 3 (blue curve), 
increases significantly, attaining a magnitude comparable 
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FIG. 22: (Color online). Persistent current at a particular AB 
flux (0 — 0.25) as a function Rashba SO interaction strength 
for a disordered {W = 2) half-filled ring with N = 60 when 
toso is set to zero. 

to that in a perfectly ordered ring. It is to be noted that 
the strength of the RSOI is strongly dependent on gate 
voltage. An enhancement of the persistent current in 
the presence of disorder can be achieved even with much 
lower values of the RSOI parameter compared to what 
have already been presented in the figures. To achieve 
this one needs to increase the size of the mesoscopic ring. 
We have checked this with a 100-site ring where even 
with tjiso = 0.5 the current increases by an order of mag- 
nitude compared to the case when t/jso ~ 0. However, we 
present the results using a somewhat larger values of tj^so 



for a better viewing of the results. Similar observations 
are made by setting tjiso = and varying t^so. 

Disorder introduces quantum interference which leads 
to localization of the electronic states. RSOI, on the 
other hand, introduces spin flip scattering in the sys- 
tem, which can destroy quantum interference effect, lead- 
ing to a possible delocalization of the electronic states. 
This leads to an enhancement of the persistent current 
in the presence of disorder. The competition between 
the strength of disorder and the RSOI is also apparent in 
Fig.[22l For small values of the RSOI, the disorder domi- 
nates. As the strength of the RSOI is increased, the spin 
flip scattering starts dominating over the quantum inter- 
ference effect, and finally the oscillations become quite 
similar to that in a ballistic ring. As the SO interaction 
is a natural interaction for a quantum ring grafted at 
a heterojunction, we are thus tempted to propose that 
the spin-orbit interaction is responsible for an enhanced 
persistent current in such mesoscopic disordered rings. 

Before we end this section, we would like to mention 
that the presence of DSOI alone leads to exactly similar 
results as expected, since the Rashba and the Dressel- 
haus Hamiltonians are related by a unitary transforma- 
tion. This does not change the physics. We also examine 
the behavior of persistent current in presence of both the 
interactions. The amplitude of the current does not in- 
crease significantly compared to the case where only one 
interaction is present. However, the precise magnitude 
of the current is sensitive to the strength of the magnetic 
flux threading the ring. The observation remains valid 
even when the strengths of the RSOI and DSOI arc the 
same. 



V. SUMMARY AND CONCLUSIONS 

In this review we have demonstrated the magneto- 
transport properties in single-channel rings and multi- 
channel cylinders based on the tight-binding framework. 
Several anomalies between the theoretical and experi- 
mental results have been pointed out and we have tried 
to remove some of these discrepancies. The main con- 
troversy is associated with the proper determination of 
persistent current amplitude. We have addressed two 
different possibilities of getting enhanced current ampli- 
tudes and proved that the currents are quite comparable 
to the experimentally predicted results. In one approach 
we have concluded after an exhaustive numerical calcu- 
lation that the presence of higher order hopping integrals 
and electron-electron interaction can provide a persistent 
current which may be comparable to the actual measured 
values. We have justified these results both for the single- 
channel and multi-channel cases. In other approach we 
have established that in presence of spin-orbit interaction 
a considerable enhancement of current amplitude takes 
place, even for the non-interacting nearest-neighbor hop- 
ing model, and the magnitude of the current in a dis- 
ordered ring becomes almost comparable to that of an 
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ordered one. In addition to these, we have also analyzed 
the detailed band structures, low-field magnetic suscep- 
tibilities, variation of electronic mobility and some other 
issues to get the full picture of the phenomena at the 
meso-scale and nano-scale levels. 

Here we have considered several important approxi- 
mations by ignoring the effects of temperature, electron- 
phonon interaction, etc. Due to these factors, any scat- 
tering process that appears in the systems would have 
influence on electronic phases. At the end, we would like 
to say that we need further study in such systems by 
incorporating all these effects. 

Future directions and opportunities: Although the 
studies involving the mcsoscopic rings and cylinders have 



already generated a wealth of literature there is still need 
to look deeper into the problems both from the point of 
view of fundamental physics and to resolve a few issues 
that have not yet been answered in an uncontroversial 
manner. For example, it may be interesting to study the 
combined effect of electron-electron interaction and spin- 
orbit interaction on magneto-transport properties. Spe- 
cially, the effect of SO interaction on electron transport 
in mesoscopic and nano-scale semiconductor structures 
should be carefully examined. The principal reason is its 
potential application in spintronics, where the possibility 
of manipulating and controlling the spin of the electron 
rather than its charge, plays the all important roleS^"—. 
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